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The dark energy model with the equation of state pd ~ ~Pd — Ap^ is studied. The model 
comprises and provides realization of several types of singularities in different parameter regimes: 
the divergence of the dark energy density and pressure at finite time and finite value of the scale 
factor, the singularity of the "big rip" type and the sudden future singularity recently introduced 
by Barrow. For parameter choices which lead to a nonsingular expansion of the universe, various 
types of the asymptotic evolution are found. The entire time evolution of the universe is described 
both analytically and numerically. The advantages of this dark energy EOS as a parametrization of 
dark energy are discussed. 
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The diverse cosmological observations, including the 
supernovae of the type la (SNIa) Q^, the anisotropies of 
cosmic microwave background radiation (CMB) [J, and 
the large scale structure (LSS) 0, among others, indicate 
that the universe is presently in the state of accelerated 
expansion. Unraveling the source of this acceleration has 
become one of the highest priorities of theoretical cos- 
mology. An interesting hypothesis is the assumption of 
modification of the gravitational interaction at cosmolog- 
ical scales However, a large majority of cosmological 
models explain the acceleration of the universe in terms 
of a component with the negative pressure, the so called 
dark energy 5]. 

The nature of dark energy has not yet been established. 
Apart from the negativity of its pressure, very few of its 
properties have been unambiguously fixed by the obser- 
vational data. It is customary to model dark energy as an 
ideal homogeneous fluid with an equation of state (EOS) 



Pd = wpd , 



(1) 



where pd and pd represent the energy density and the 
pressure of dark energy, respectively. The parameter of 
the equation of state w is generally a function of cosmo- 
logically relevant quantities. No matter if the parameter 
of the EOS is considered constant or is modeled as the 
dynamic quantity, the analyses of the cosmological data 
constrain the present value of the parameter of the dark 
energy EOS to an interval around the value w = — 1 0. 

A natural benchmark for the dark energy EOS pa- 
rameter w is the value characteristic of the cosmological 
constant (CC), i.e., w = —1. This benchmark value is 
very useful for natural classification of dark energy mod- 
els. Namely, the dark energy density in the models with 
w > —1 decreases with the expansion of the universe. 
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For the models with w — —1 (CC), the dark energy den- 
sity is constant, whereas it grows with the expansion of 
the universe for the models with w < —1. The cosmo- 
logical constant is one of the most popular explanations 
of the nature of dark energy Jl . A huge discrepancy be- 
tween the value inferred from the observations and the 
value predicted by quantum field theory requires the ex- 
istence of some unknown fine-tuning mechanisms. The 
coincidence of the CC energy density and the energy 
density of nonrelativistic matter at the present epoch 
of the expansion of the universe requires further fine- 
tuning. Some of the problems related to the CC may be 
alleviated in the formalism of the renormalization group 
equation for the CC H H E IH El- This approach 
is also promising in the treatment of some astrophysical 
problems IQi, il^j . A large class of dark energy models 
satisfy w > —1. The prominent members of this class are 
quintessence k-essence 0|, and Chaplygin gas fl^ 
which incorporates a very attractive possibility of unify- 
ing dark energy and dark matter. The models satisfying 
w < —1 are called phantom energy models [T^ llTl lT8| . 
Some phantom energy models include a very interesting 
"big rip" event with a dramatic influence on the destiny 
of all bound structures 19] . Certain attempts of the mi- 
croscopic formulation of phantom energy models exhibit 
the quantum instability of vacuum '20'| , which further in- 
creases the controversy around the concept of phantom 
energy. Some recent considerations also study the possi- 
bility that the expansional effects, usually attributed to 
phantom energy, might be mimicked by the more com- 
plex dynamics of nonphantom components or configu- 
rations [231 . Finally, a legitimate question is whether 
the parameter of the EOS may cross the CC bound- 
ary w = — 1. A recent analysis '2^ shows that in the 
framework of k-essence dark energy, such a transition is 
unlikely. It should also be mentioned that interesting al- 
ternatives explaining superaccelerated expansion of the 
universe in terms of quantum effects of scalar fields in 
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the de Sitter space, requiring no phantom energy, also 
exist m. 

All intervals for w obtained as the constraints from the 
observational data include the w ~ —1 value and signifi- 
cantly overlap with both w < —1 and w > —1 segments 
of the w axis. From the theoretical side, at w = —1 
there is a qualitative change in the evolution, as well 
as in other properties of dark energy. Therefore, since 
from both theoretical and observational perspective the 
dark energy EOS is "centered" at the vacuum EOS, it 
is useful and interesting to consider dark energy mod- 
els with the EOS which is also centered at the vacuum 
EOS w = —1. In this paper we consider a class of such 
models where the dark energy EOS is reminiscent of an 
expansion around the vacuum EOS. Namel y, w e consider 
dark energy models (already mentioned in (2^) with the 
following equation of state: 



Pd = -pd - Ap2 . 



(2) 



Here A and a are real parameters. In our considerations 
we focus on the models with a ^ 1 because in the a = 1 
case the dark energy EOS reduces to the dark energy with 
the constant EOS of either the quintessence {A < 0) or 
the phantom {A > 0) type. 

Many interesting dark energy models are motivated 
by the considerations from the string theory [2^, one 
specific example being the tachyonic dark energy [26j |. 
Consideration of the contribution of the transplanckian 
degrees of freedom in the framework of string theory |23| 
implies the possibility of the supernegative equation of 
state for the dark energy and also further motivates the 
analysis of the more general dark energy equation of state 
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From the conservation of the energy-momentum tensor 
of dark energy with the EOS 10 



da 

dpd + 3(/9d +Pd) — = 



(3) 



we obtain the following expression for the scaling of the 
dark energy density: 



Pd = Pdfi 1 + 3ii(l - a) In ■ 



ao 



1/(1-") 



(4) 



Here we have introduced an abbreviation A = Ap^^^. 

The parameter A determines the type of the behavior of 
the dark energy density with the expansion of the uni- 
verse. The speed of change of the dark energy density 



J / \ i/(i-")-i -1 

^ = 3ip,,o l + 3i(l-a)ln-'' ^ 
da \ ao 



(5) 



shows that for A > the dark energy density grows, for 
A = the dark energy is constant, while for A < it 
decreases with the expansion of the universe. 

We assume that the universe contains two components, 
the dark energy with the EOS Q and the nonrelativistic 



matter. The evolution of the universe is determined by 
the Friedmann equation 



{pd + Prn): 



(6) 



where dots denote time derivatives and pm stands for 
nonrelativistic matter energy density. 

In the analysis of the dynamics of the universe we first 
focus on the fate of the universe, i.e., its asymptotic ex- 
pansion. We assume that the dark energy density in- 
creases, remains constant or decreases more slowly than 
the nonrelativistic matter energy density or the curva- 
ture term fc/a^. The time dependence of the scale factor 
for a 7^ 1/2 is then given by the expression 



1-2q 

a, \ 2(i-Q) 
1 + 3A{1 - a)ln — 
ao, 

l-2a 
2(1- = ) 



l-f^3i(l-a) In — 
|i(l - 2a)nl(^Ho{h - h) 



(7) 



Here fld,o = Pd,o/Pc.,o, where pc.o = 3iJo/(87rG') is the 
critical energy density at the present epoch. 

For a = 1/2, the expression for the dependence of the 
scale factor on time becomes 



In- 



1 



|Aln^ 

2 ao 



1 



t2). 



(8) 



Once we have the complete set of equations determin- 
ing the temporal evolution of the universe, we can study 
the parametrical dependence of the cosmological solu- 
tions. Equations Q), 0, and © single out two specific 
values for the parameter a: a ~ 1/2 and a = I. These 
values represent the bordering points of the intervals on 
the a line with qualitatively the same cosmological be- 
havior. Therefore, in our analysis of the parameter space 
of the model, for each value of the parameter A we con- 
sider the following four intervals for the parameter a: 
a > 1, 1 > a > 1/2, a = 1/2, and a < 1/2. The value 
a = 1 is not analyzed since it corresponds to the famil- 
iar model of dark energy with the constant parameter 
of the EOS. We first consider the positive values of the 
parameter A which correspond to phantom energy. 

A > 0, a > 1. Equation Q shows that, with these 
parameter values, the pressure of dark energy grows much 
faster compared with phantom energy with the constant 
parameter of the EOS. On the basis of this observation, 
one might expect even more singular behavior than in 
the usual "big rip" case. The expression for the scahng 
of the dark energy density with the scale factor a, more 
conveniently written as 



Pd = Pd,o 1 - 3A{a - 1) In — 
* ao 



-l/(a-l) 



(9) 
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FIG. 1: The time evolution of the scale factor of the universe 
for ild,o = 0.7, Qm,o = 0.3, A = 1 and four typical values of 
the parameter a: a = —1,0.5,0.75,2. For a = 2, the scale 
factor reaches the finite value at the singularity while for a — 
0.75, the scale factor of the universe diverges at singularity. 



FIG. 2: The time dependence of the Hubble parameter for 
A^l, Qdfi = 0.7, Qmfi = 0.3, and a = -1,0.5,0.75,2. For 
a = 0.75 and a = 2 there is the singularity of the Hubble 
parameter appearing at finite time. 



shows that the dark energy density diverges for the finite 
value of the scale factor 



floe 



l/(3A(a-l)) 



(10) 



At the same value of the scale factor the dark energy 
pressure diverges towards negative infinity. Equation (jj)) 
reveals that this singularity is reached at finite time tmax- 
The evolution of the scale factor before the singularity 
can be expressed as 



aoexp 



3A{a - I 



(1- 



lAi2a- 



t) 



2q-1 



(11) 



To summarize, the evolution of the universe reaches the 
singularity in which the dark energy density and the pres- 
sure diverge at the finite time tmax and the finite value 
of the scale factor amax- 



A > 0, 1/2 < a < 1. Equation igj shows that in this 
case the model exhibits the minimal value of the scale 
factor 



(12) 



Let us denote the instant of time when the description 
of the evolution of the universe in terms of dark energy 
only becomes appropriate by t.^, and the corresponding 
scale factor by a*. Equation {Tj) acquires the form 



l + 3A(l-a)ln — 
ao 



'2(l-c) 



3 



-^Ai2a-l)nll^Hoit~U) 



2(1-^) 
2q-1 



(13) 



At the instant of time 

_ 2a- 

2 ( ~ a \ ^'1' 

t„p = U + — ^ 1 + 3.1(1 - a) In — 

3A{2a-l)nY^Ho \ J 

(14) 

the scale factor diverges. The dark energy density and 
pressure diverge at the same instant. This scenario is 
qualitatively equivalent to the "big rip" scenario present 
in many phantom energy models. 



A > Q, a = 1/2. Equation JSJ gives the law of evolu- 
tion 



a = aoexp 



2 



3 , a* 
1 + - Aln — I X 

2 ao. 



exp[^Anl(^Ho{t-U)] -1 



(15) 



Here the quantities and a* have the same inter- 
pretation as in the preceding case. The model ex- 
hibits the minimal value of the scale factor amm = 
aoexp{—2/{'3A)). The future expansion of the universe 
has no singularities and continues infinitely. 

A > 0, a < 1/2. In this case, the model also has a 
minimal value of the scale factor. The future evolution 
of the universe has no singularities and is given by the 
expression 
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FIG. 3; The time evolution of the dark energy density for 
A = 1, fld.o ~ 0.7, Qm.o = 0.3, and the typical values a = 
— 1, 0.5, 0.75, 2. For a = 0.75 and a = 2, there is the onset of 
singularity of the dark energy density at finite time. 
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FIG. 4: The time dependence of the dark energy pressure 
for A = 1, fld,o = 0.7, Qm.o = 0.3, and the typical values 
a = —1, 0.5, 0.75, 2. For a = 0.75 and a = 2, the dark energy 
pressure diverges to — oo. 



(16) 



The results on the asymptotic expansion of the uni- 
verse, obtained in an analytically tractable form, can be 
further supported and completed by numerical calcula- 
tions. The full time dependence of the universe expansion 
obtained numerically is depicted in Figs HI El and ^ 

The analysis presented so far has been concentrated on 
the models with positive A which are of phantom char- 
acter. Next we focus on the models with negative A. In 
these models the character of dark energy is nonphantom. 



The accelerated expansion of the universe at the 
present epoch implies the negativity of the dark energy 
pressure at present. This requirement was automatically 
satisfied for the models with the phantom character. For 
models with A < 0, this requirement produces nontriv- 
ial constraints on the parameter A. Namely, if we re- 
quire the negativity of the dark energy pressure at the 
pressent epoch, the parameter A must obey the con- 
straint A > —1. If we demand the negativity of the dark 
energy pressure at some z > 0, an analogous condition 
for A can be obtained in a straightforward manner and it 
is generally a dependent. The behavior of nonphantom 
dark energy models is further classified with respect to 
the parameter a. 

A < 0, a > 1. The scaling of the dark energy density 
with the scale factor 

/ , ^x-r/(a-r) 
Pd = Pd,o (l + 3A(l-a)ln— j , (17) 

points to the existence of the minimal value of the scale 
factor amin = aoexp{—l/{3A{l — a))). The energy den- 
sity pd decreases with the expansion of the universe. 



The term proportional to p" becomes less significant 
with the expansion of the universe and the dark energy 
EOS asymptotically approaches the cosmological con- 
stant EOS. These findings on the scaling properties of 
the dark energy density are also refiected in the asymp- 
totic evolution of the scale factor: 



aoexp 



1 



3A{1 - a) 



1 + 3A{1 - a)ln- 



^i(l ~ 2a)nl/^Ho{t ~ U 



1-2q 



(18) 



The future expansion of the universe is unbounded and 
nonsingular. 



We continue the discussion with some general remarks 
for all cases with a < 1 since they share the same type 
of behavior of the dark energy density. Namely, the dark 
energy density decreases and at the scale factor 



ciNULL = floe 



-l/{3A{l-a)} 



(19) 



the dark energy density vanishes. For these parameter 
values, the expressions for the asymptotic expansion are 
no longer applicable since dark energy does not domi- 
nate the future expansion of the universe. From the dark 
energy EOS Q we can see that at a^uLL the dark en- 
ergy pressure vanishes for < a < 1, acquires the finite 
value pd = —A for a = and diverges for a < 0. Let us 
describe the dynamics of the universe in more detail for 
different intervals of the parameter a. 

A 1/2 < a < 1. Let us discuss the future behavior 
of the universe for a simplified cosmological model with 
a single component which is dark energy with the EOS 
In this setting, equations ([TJl and ijHJl become exact 
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FIG. 5: The time evolution of the scale factor of the universe 
for fld,o = 0.7, flm,o = 0.3, A — —0.5 and a = —1. The scale 
factor of the universe is finite at the onset of singularity. 



for the flat universe. Although the value unull is finite, 
it is reached only at infinite time, as can be seen from l(7|l. 
Therefore, in the model with the dark energy component 
only, there is no finite time singularity. The quantities 
Pd and Pd asymptotically go to 0. 

A < 0, a = 1/2. In the model with the dark energy 
component only, introduced in the preceding paragraph, 
equation ((SJ governs the expansion of the universe. The 
expansion reaches aNULL at infinite time and there is no 
finite time singularity. Both pd and pd asymptotically 
vanish. 

A < 0, < a < 1/2. In the model with the dark 
energy component only, the universe reaches the finite 
scale value unull at finite time tNULL = —2/{3A{l — 

1/2 

2a)Qj Hq), as follows from iQ. At Inull the quantities 
Pd and Pd vanish. The expansion of the universe stops 
with a(tjv[/LL) = and a{tiqijLL) = 0. There is no 
singularity at the finite time ^jvc/ll- 



In the interval < a < 1, the scaling of pd with 
a is generally not well defined by (QJ for a > aj^iuLL. 
However, for countably many a values the expression 
Q can be extended beyond a null- Namely, for val- 
ues a = 1 — l/(2n), where n — 1,2,..., the dark energy 
density is well defined for all scale factor values. If we al- 
low the possibility of negative values for the dark energy 
density, for a = 1 — l/{2n -1-1), where n = 0, 1, 2, . . ., the 
dark energy density is well defined for all scale values as 
well. However, for negative values of p^, the dark energy 
pressure is not well defined. 

The four preceding paragraphs discuss the evolution of 
the universe in a model containing only the dark energy 
component with the EOS (0). A realistic model with 
both nonrelativistic and dark energy components can be 
easily numerically handled for scale factors up to unull- 
Furthermore, the numerical analyses show that unull is 



reached at finite time for parameter values < a < 1. 
At unull the quantities pd andpd vanish, but pm and H 
do not. At a^uLL there is no singularity, which is an ex- 
pected feature. Namely, in the situation where the dark 
energy component does not cause any singularity, it is 
highly improbable that the addition of the well-behaved 
nonrelativistic component would induce singular behav- 
ior. The situation at a^uLL implies that the expansion 
should continue beyond a^uLL- However, at a > qnull 
the scaling of pd with a is not well defined Q). It is 
important to notice that the scaling equation Q is also 
satisfied with pd = (and pd = 0). Therefore, a possible 
explanation of the expansion at a > unull is that once 
the dark energy density reaches at a^uLL, it remains 
in the state with = and pd = Q. The expansion 
at a > ttNULL is entirely governed by the nonrelativis- 
tic matter and the dark energy is gravitationally ineffec- 
tive. This interpretation yields a cosmological model in 
which the expansion at a > ajsnjLL is decelerated, and 
(for a flat universe) the expansion continues to infinity. 
Furthermore, during the entire expansion, there are no 
singularities in a and a. 

^ < 0, a < 0. The dark energy density tends to 
= as the scale factor approaches a^uLL- At the 
same time the pressure diverges, i.e., pd ^ oo. In 
the model which contains only the dark energy com- 
ponent, this singularity is reached at a finite moment 
tNULL = -2/(3i(l - 2a)D}-J'^Ho). The numerical cal- 
culations for the more realistic model with both nonrel- 
ativistic matter and dark energy components also show 
that the singularity is reached at finite time. The time 
dependences of relevant cosmological quantities arc dis- 
played in Figs El and |S1 The characteristics of this 
singularity are as follows: the singularity is encountered 
at finite time; at the singularity the scale factor is finite. 
Fig. 13 the quantities a and H are finite. Fig. El and so 
are the energy densities pd (Fig. [7|) and ptot = Pm+ pd] 
the acceleration and the dark energy pressure (which 
equals the total pressure) diverge (Fig. jSJ 

a —f — oo , Pd +00 . (20) 

The characteristics of this singularity match those of the 
sudden future singularity introduced by Barrow [2^ . 
Therefore, the cosmology with the dark energy which 
has the EOS @ for A < Q and a < leads to the fi- 
nite sudden singularity which satisfies the strong energy 
condition. This example shows that a relatively simple 
dark energy EOS may lead to the type of singularity pro- 
posed by Barrow. Some other types of cosmological sin- 
gularities characterized by an infinite deceleration have 
recently been introduced 30]. 

From the overview of the different types of the expan- 
sion of the universe, we can see that in various parameter 
regimes the dark energy model Q leads to various types 
of singularities or asymptotic behaviors. The wealth of 
phenomena embodied in one model makes it an attractive 
object of theoretical study. At this moment we would like 
to emphasize that the model has some advantages beyond 
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FIG. 6: The dependence of the Hubble parameter on time for 
fld,o = 0.7, Qmfi = 0.3, A = -0.5 and a = -1. The Hubble 
parameter is finite at the onset of singularity. 



FIG. 7: The dark energy density dependence on time for 
fld,o = 0.7, n™,o = 0.3, A = -0.5 and a = -1. The dark 
energy density vanishes at the singularity. 



theoretical considerations. Namely, dark energy with the 
EOS (O is very suitable for the analysis of cosmologi- 
cal data. The dark energy EOS is centered around the 
cosmological constant EOS and the parameters A and 
a measure the deviation from the cosmological constant 
EOS. The sign of the parameter A determines vifhether 
dark energy is of the phantom or the quintessence char- 
acter. Furthermore, the parameter A determines the am- 
plitude of deviation from the cosmological constant EOS, 
while the parameter a determines the speed of variation 
of this deviation. It could be said that the parameters A 
and a measure the deviation from the ACDM cosmology. 
Different intervals of the parameter a, for the fixed A, 
imply the qualitatively different destiny of the universe. 
The constraint from the observational data to any of the 
intervals for the parameter a could qualitatively deter- 
mine the fate of the universe and eliminate some types of 
singularities or asymptotic behaviors. Finally, this dark 
energy model provides a natural framework for the time- 
varying parameter of the EOS. All these characteristics 
put forward the dark energy model (0) as a promising 
candidate for the analysis of cosmological data. 



The analysis presented so far is based upon the ex- 
pression ^ as an exact EOS. However, the usefulness 
of the results of this analysis surpasses the model itself. 
It would certainly be of interest to investigate the mod- 
els obtained by adding new terms into the dark energy 
EOS (01. This would lead to a general dark energy EOS. 
Even in a general dark energy model the results of this 
paper are still useful, since the EOS lO can in this case 
be considered as an approximation of the more general 
EOS around the present epoch of the evolution of the 
universe. In this case, the solutions obtained above de- 
scribe the expansion of the universe around the present 
epoch well, whereas the description of the singularities 



obtained in the analysis of Q is generally no longer ap- 
plicable. The analysis of the dark energy model presented 
in this paper reveals that even a relatively simple dark 
energy EOS may comprise very different types of singu- 
larities. It is reasonable to expect that for a model with 
a more general EOS we may expect even much richer 
abundance of phenomena. A very concrete contribution 
of the EOS (0) is the demonstration which EOS may lead 
to the type of singularity at finite time and finite scale 
factor and to an opposite type of the sudden future sin- 
gularity |28j |. Another important conclusion is that both 
of these singularities follow from the same dark energy 
EOS, but in different parameter regimes. 



The considerations of the model with the EOS (|2Jl pre- 
sented so far have been devoted to the future evolution 
of the universe. Let us briefly discuss the consequences 
of the dark energy model considered in this paper for 
the past evolution of the universe and the constraints on 
the model parameters which follow from the comparison 
againts the data obtained from various cosmological ob- 
servations. For A > and a > 1 the dark energy density 
grows with the expansion of the universe, as does the 
absolute value of the dark energy pressure. Therefore, 
at large redshifts the dark energy density and pressure 
become negligible compared to the energy densities of 
other components of the universe and negligibly affect 
the physics of the early universe (such as primordial nu- 
cleosynthesis). For the case ^ > and a < 1 there is the 
minimal value of the scale factor amin (|12|l at which the 
dark energy vanisehes and is not well defined for smaller 
scale factors. Furthermore, for a < 0, the absolute value 
of the dark energy pressure becomes very large and di- 
verges at amin- In order to avoid significant impact on 
the early universe dynamics, the minimal scale factor of 
the model needs to be very small. From H12|) follows 
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FIG. 8: The dark energy pressure dependence on time for 
fld,o = 0.7, Qmfi = 0.3, A = -0.5 and a = -1. The dark 
energy pressure diverges towards +oo. 



that A{a — 1) must be a very small number. In the case 
A < and a > 1 there is the minimal value of the scale 
factor a„jj„ again. In this case, however, the dark en- 
ergy density diverges. Clearly, amm has to be very small 
in order to maintain the consistency with the results of 
the early universe physics. This again restricts the prod- 
uct A{a — 1) to be a sufficiently small number. For the 
case A < and a < 1, the dark energy density diverges 
as the scale factor goes to and since its dependence 
on the scale factor is of the logarythmic type, no severe 
constraints are expected for finite scale factor values. Fi- 
nally, it is important to say that these considerations on 
the past behavior of the universe are made under the 
assumption that the the dark energy EOS is valid even 
for the smallest values of the scale factor. If one con- 
siders ((SJ only as a late time effective dark energy EOS, 
the conditions on parameters a and A may be much less 
stringent. 



The treatment of this paper is entirely classical. In 
the vicinity of singularities the quantum effects come into 

t, i.e., the universe enters a new quantum gravity era 
It has been shown that quantum effects may 
reduce or even completely remove the singular behavior 
at classical singularities. In that respect, the study of 
the interplay of classical and quantum effects is needed 
to obtain a complete physical picture of the dynamics 
near the classical singularity. 

In conclusion, the dark energy model with the EOS 
exhibits a number of interesting phenomena in dif- 
ferent parameter regimes. The model comprises several 
types of singularities appearing at some finite time after 
the present epoch of the evolution of the universe. For 
A > and a > 1, the universe encounters the singu- 
larity at finite time and finite value of the scale factor, 
where the dark energy density and pressure diverge. For 



A > and 1/2<q;<1, at singularity, which is reached 
at finite time, the scale factor and the dark energy den- 
sity and pressure diverge. This singularity is in these 
respects equivalent to the "big rip" singularity found for 
phantom energy models with the constant parameter of 
EOS. For A < and a < 0, the singularity is reached 
at finite time and finite value of the scale factor. At 
the singularity the dark energy density vanishes and the 
total energy density of the universe is finite. The dark 
energy pressure diverges to -l-oo at the singularity. This 
form of singular behavior is of the type of sudden future 
singularity introduced by Barrow 28]. The model com- 
prises a number of singularity types with very different 
properties. Furthermore, it provides a relatively simple 
framework for the concrete realization of the sudden fu- 
ture singularities such as those introduced by Barrow. In 
the rest of parameter space, the expansion of the uni- 
verse is nonsingular with, however, different asymptotic 
behavior. The model itself is a useful testing ground 
for the study of singularities in dark energy cosmologies. 
The dark energy with the EOS Q has an additional ad- 
vantage that the restrictions on the parameters A and 
a obtained by comparison with the observational data 
may qualitatively select the destiny of the universe. The 
model also provides a natural framework for the study 
of dark energy with a time-dependent parameter of the 
EOS. The features of the dark energy model Q imply 
that for some more general dark energy EOS one might 
expect an even larger abundance of interesting phenom- 
ena and singularities discussed in the literature may ap- 
pear for the relatively simple dark energy EOS. The com- 
parison with the observational data will hopefully deter- 
mine which destiny of our universe is the most reasonable 
to expect. 

Note added. In a recent paper appearing shortly 
after our paper, a general analysis of the types of singu- 
larities in phantom energy models was made. The singu- 
larity at the finite time and the finite value of the scale 
factor, first introduced in our paper, was classified as the 
type III singularity in |3lj |. Furthermore, using the inte- 
grated conformal anomaly, an analysis of quantum effects 
near singularity, was performed in [3 If . It was shown that 
for the singularity of the type III appearing for the model 
(PI the behaviour near the classical singularity becomes 
less singular. In particular, the quantum effects near the 
type III singularity may completely remove the singu- 
larity in the behavior of the Hubble parameter H. The 
reduction of the singular behavior owing to the quantum 
effects was also found for other types of singularities in 
phantom energy models. 
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